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Results of our next-to-leading order analytic calculation for longitudinally polarized photopro- 
duction of heavy quarks are presented. Certain differential cross sections and the corresponding 
asymmetries are computed for energy ranges accessible to CERN and HERA. We argue that 
determination of the polarized gluon distribution is possible at both places: at HERA certain 
\Q ', acceptance cuts are necessary. 

o 

Deep Inelastic Scattering (DIS) of longitudinally polarized particles has been very useful 
in extracting the polarized valence distributions. However, it is quite limited in extracting 
the polarized gluon distribution Ag. Much was already said at this Workshop about the 
importance of knowing Ag, and there is no need to talk about it any further. In order to 
extract Ag with a good accuracy one needs experiments involving longitudinally polarized 
o ■ reactions dominated by subprocesses with initial gluons. 

One such reaction is the inclusive polarized photoproduction of heavy quarks: 



7 + P^Q(Q)+X, (1) 

which at leading order (LO) proceeds through the subprocess: 

1 + g^Q + Q (2) 

There are several experimental proposals on the subject in various stages of approval M. 

Next-to- leading order corrections (NLOC) are very important for several reasons: a) 
The leading order estimates depend strongly on the a priory unknown renormalization and 
factorization scales; in general the NLOC reduce this uncertainty. More precisely, going 
to one level higher in the strong coupling a s may reduce the uncertainty by order a s ; b) 
NLOC may be large and negative in some regions of phase space, making cross sections 
small and difficult to measure; c) At NLO, subprocesses unrelated to the LO ones could 
come into play. In principle, the contribution of such subprocesses could be significant in 
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certain important parts of phase space and this could affect even the shape of the final 
cross sections. However, admittedly in the present case, the subprocesses 

i + q(q)^Q + Q + q(q), (3) 

(q =light quark) which belong to this class, make a relatively small contribution. The 
same holds for the corresponding unpolarized subprocesses f2|, as well as for other such 
subprocesses (taken separately) we have studied [[|. 

Nevertheless, irrespective of relative sizes, this talk presents a complete NLO calculation 
of polarized photoproduction of heavy quarks. Some results, in a very preliminary stage, 
have already been reported [|J. 

An independent calculation of NLOC has appeared and we subsequently comment 
on it. 

Let Mj(Ai, A 2 ) be the amplitude of any of the contributing subprocesses and Aj, A 2 the 
helicities of the initial partons. Then unpolarized and polarized, cross sections correspond 
to the quantities: 

±£[M i (++) Mj (++) ± M t (+-) Mj (+-)], (4) 

where S denotes summation over the helicities and colors of the final particles and average 
over the colors of the initial. Note for the unpolarized case one usually averages over n — 2 
spin degrees of freedom for every incoming boson, which amounts to multiplying (||) by 
2/(n — 2) for each initial boson. 

Note that the Abelian part of the NLOC to our dominant subprocess (fj) is of interest 
in itself, for it determines the NLOC to 7 + 7 — > Q{Q) + X ||, [5J. Note also that, as it was 
already done in ||, a calculation of NLOC for a polarized reaction and the corresponding 
unpolarized one can be done by simply using different projection operators. 

Our regularization scheme is dimensional reduction (RD) properly supplemented by 
finite counterterms || and conversion terms || (see also below). We work in the Feynman 
gauge. 

In determining the loop contributions, the renormalization of the heavy quark mass 
and wave function were carried on shell. Furthermore, in our scheme, the contribution 
of a heavy quark loop in the gluon self-energy is subtracted out, i.e. the heavy quark is 
decoupled. This is consistent with parton distributions of which the evolution is determined 
from split functions involving only light quarks. 

Our loop contributions have been determined by REDUCE || and to some extent by 
FORM PH|| . Tensor integrals have been reduced to scalar ones. It is worth noting that, in 



general, the coefficients of pole terms 1/e 2 and 1/e for separate graphs are not proportional 
to the LO term; only their sums are. 

Gluon Bremsstrahlung contributions of e.g. 7 + p — > Q + X are determined, as usual, 
by working in the cm. frame of Q and g (Gottfried- Jackson). Many necessary n- and 
4-dimensional integrals are given in [Tl[] and the rest were derived by ourselves ||. Fac- 



torization terms are added in the usual way, and the cancellation of the loop 1/e 2 and 1/e 
singularities is a good indication of the correctness of our results. 
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Figure 1: Differential asymmetry as defined in the text for three sets of GS polarized pd. Solid 
lines: set A, dashed: set B, dotted: set C. 



Here we leave out a resolved 7 component; this, as well as an estimate of the anticipated 
experimental errors, is given in ||. 

Our cross sections are convoluted with parton distributions evolved via two-loop split 
functions. However, polarized split functions have been determined in the modified t'Hooft- 
Veltman (HV) scheme [ 12fl . Thus, conversion terms were necessary to transform our results 



to the HV scheme. 

We now present some numerical results for the polarized differential cross sections: 
Ada/dp? and Ada/dy, where px an d y the transverse momentum and cm. rapidity with 
respect to the photon of a heavy quark of mass m. For the renormalization (/i) and 
factorization (M) scales we take: /i 2 = M 2 = p\ + m 2 for Ada/dpx and /1 2 = M 2 = 4 m 2 
for Ada/dy. With y/S the cm. energy of the colliding particles we define xt = 2px/\fS- 
We use the three sets of polarized parton distributions (pd) of ||13|| , which mainly differ in 
Ag: Sets A and B differ in the size of Ag, set C also in the shape. Our asymmetries are 
defined as follows: 



For the unpolarized cross sections we use the pd of [14]]. Our figures are for three energies: 



VS = 10 GeV (COMPASS and SLAC) and = 100 and 300 GeV (HERA). 

Fig. 1 presents our NLO All (pt) for Q = c and Q = b quark. Clearly, in general, sets 
A, B and C lead to Ada/dp^ s well separated. Also note that, in general, at small xt 



3 




-1.5 -1 -0.5 0.5 1 1.5 



y 



Figure 2: Rapidity differential asymmetry as defined in the text for three sets of GS polarized 
pd. Lines as in Fig. 1. 

s A LL (p T ) are relatively small; exception is Vo = 10 GeV, when A LL (p T ) stay relatively 
sizable. On the other hand, at low xt both the polarized and unpolarized cross sections 
are relatively large; this is particularly true for the unpolarized ones at HERA. Then, at 
HERA, at small xt, A ll (p T ) are very small, as was indeed found in ||. 

Fig. 2 presents All (y) and Fig. 3 the sizes of Ada/dpx and Ada/dy. In the later we 
do not plot results for v^S" = 300 GeV as the cross sections are too small to be measured. 
For c-quark it seems that Ada/dpx is measurable up to xt — 0.4. Also, at v^S* = 10 GeV, 
looking at Fig. 2 and Fig. 3 we see that in the range 1 < y < 1.5 both A LL (y) and Ada/dy 
are quite sizeable, so useful information on Ag could well be obtained. 

Ada/dy of Fig. 3 show also that heavy quarks are mainly produced forwards with 
respect to the (incoming) photon. 

As we discuss in detail in ||, in general, we agree with M. Nevertheless, perhaps on 
the basis of the smallness of HERA All at small xt, Ref. concludes that HERA is 
rather useless in specifying Ag. Here we argue that it may be not so: On the basis of 
Fig. 1, reconstruct events and select only those with, say, xt > 0.2; in other words, carry 
integrations of Ada/dpx over some cut phase space. This may well enhance the resulting 
All- 

Of course, to reach a definitive conclusion, more work is needed, including an estimate 
of the corresponding errors. 

We would like to thank V. Spanos and G. Veropoulos for participating in part of the calcu- 
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Figure 3: Polarized cross sections for two sets of GS polarized pd. Solid lines: set A, dashed: set 
B, dotted: contribution due to the light quark subprocess (||). Leading order results are marked 
with stars. 
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